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Nuclear site analysis methods are used to enumerate the normal modes of ABX3 perovskite poly- 
morphs with octahedral rotations. We provide the modes of the fourteen subgroups of the cubic 
aristotype describing the Glazer octahedral tilt patterns, which are obtained from rigid rotations of 
the BXq octahedra with different sense and amplitude about high symmetry axes. We tabulate all 
normal modes of each tilt system and specify the contribution of each atomic species to the mode 
displacement pattern, elucidating the physical meaning of the symmetry unique modes. We have 
systematically generated 705 schematic atomic displacement patterns for the normal modes of all 
15 (14 rotated + 1 unrotated) Glazer tilt systems. We show through some illustrative examples 
how to use these tables to identify the octahedral rotations, symmetric breathing, and first-order 
Jahn- Teller anti-symmetric breathing distortions of the BXq octahedra, and the associated Raman 
selection rules. We anticipate that these tables and schematics will be useful in understanding the 
lattice dynamics of bulk perovskites and would serve as reference point in elucidating the atomic 
origin of a wide range of physical properties in synthetic perovskite thin films and superlattices. 
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I. INTRODUCTION 

The compounds with the general formula ABX3 (X = 
or a halogen) where A is either an (alkaline) earth or a rare 
earth metal and B is a transition metal with partly filled 
d-orbitals form the perovskite crystal family. Perovskites 
are a truly versatile material class that is malleable both 
chemically and structurally. At least 30% of the elements 
from the periodic table can occupy the A cation position 
and over half of the periodic table fills the B cation 
position, with almost 100% substitution^ The structural 
adaptability of the ABX3 perovskites, enabled by the 
flexible and corner-connected BXq octahedral network, 
explains why such diverse chemistries are compatible in 
the perovskite crystal structure (Fig. [T]) . 

The chemical and structural compatibility of oxide- 
based perovskites at the atomic-scale makes it possible to 
tune their macroscopic properties in a variety of ways, e .g., 
by judicious choice of the cations (chemical pressure) pEl 
and hydrostatic pressure.^ Consider the perovskite ox- 
ide series AN1O3, where A is a rare-earth element: The 
Ni-O-Ni bond angles between neighboring octahedra can 
be made to approach 180° (cubic symmetry), driving an 
insulator-metal transition as the A-site ionic radius is 
increased systematically by substitution of Eu, Nd, Pr 
and LaP A modern technique to tune structure-derived 
electronic properties in perovskites includes octahedral 
engineering, where substrate-imposed strain and superlat- 
tice formation are used to alter the magnitude and flavor 
of the BOq octahedral rotation patterns by exploiting 
both misfit strain-octahedral rotation and heterointerfa- 
cial coupling between two similar or dissimilar tilt pat- 
terns possessed by the substrate and filmPZI The directed 
control of the crystal structure in epitaxial films often 
stabilizes states which are not found in bulk equilibrium 
phases, producing dramatic changes in their macroscopic 



properties. 

Control over macroscopic properties requires an under- 
standing of the atomic structure and an appropriately 
facile technique to measure both the property-dictating 
structural units - in this case the BOq octahedra - and 
their response to external stimuli. In other words, appro- 
priate experimental methods are required to determine 
the structural transitions in functional perovskite oxides 
induced by hydrostatic or chemical pressures, strain and 
superlattice formation. However, the measurement of 
oxygen positions in ultrathin ABO3 films, required to 
identify the octahedral tilt patterns, is non-trivial for 
a variety of experimental reasonsP Thin films present 
limited sample volume and often are grown on relatively 
thick substrates, which in addition to weak scattering 
from the oxygen atoms, makes structure determination 
of the film challenging. Raman s pect roscopy, along with 
real and reciprocal space imagingJSl^J however, are experi- 
mental techniques that can reliably determine octahedral 
tilt patterns. 

Raman scattering by zone-center, k — (0, 0, 0), phonons 
is an established method for identifying structural phases 
in solid, thin-film, and nanostructured semiconductor 
and oxide materials, since vibrational modes provide a 
unique signature of the crystal structure. Analyses of 
the appearance, energies, relative intensities, linewidths 
and lineshapes, and symmetry of Raman peaks collected 
have shown Raman scattering to be indispensable for 
investigating important relationships among correlated 
structural, electronic charge and spin, and dipolar de- 
grees of freedom in complex oxide materials.^ Despite 
the optical diffraction-limited spatial resolution, Raman 
scattering can probe the onset of a new phase with lo- 
cal correlation lengths of a few unit cells. Furthermore, 
judicious use of the optical polarization vectors' relation- 
ship with respect to the plane of the sample, which is 
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obtained through cross sectional sample preparation, can 
selectively enhance the Raman scattering intensity by 
selected phonons, particularly in the vicinity of a phase 
phase transition. 

Raman spectroscopy studies on cubic ABO3 type per- 
ovskite systems (e.g., BaTiOs) have been carried out by 
DiDomenico et al, focusing on the the splitting of the 
Raman active modes due to the cubic-to-tetragonal trans- 
formation in these systems.^ More recently, Raman spec- 
troscopy on ABO3 type perovskite oxides have focus ed on 
the orthorhombic and the rhombohedral structures! 13 * 14 ! 
The primary displacement patterns investigated in these 
studies include (i) the soft modes associated with the 
rotation of the rigid BOq octahedra about different crys- 
tallographic axes, and (ii) the higher energy symmetric 
and anti-symmetric octahedral breathing-type modes as- 
sociated with first-order Jahn- Teller effects. 

Dubroka et al. have shown in L&i-ySr y MTii^ x M x Fe03 
(A/ = Cr, Co, Cu, Zn, Sc or Ga) that the frequency of the 
soft mode evolves linearly with the octahedral rotation 
angle about the [110]-pseudocubic direction with a pos- 
itive slope of 20 cm -1 per degree-rotation.^ Since the 
oxygen positions dictate the octahedral tilt angle, Raman 
scattering provides a non-destructive and readily acces- 
sible means to determine the oxygen positions in ABO3 
perovskites, for which, alternative experimental methods 
require substantial investments in sample preparation, 
equipment and/or data processing time. 

Higher-frequency modes - the Jahn- Teller activated an- 
tisymmetric breathing modes and the symmetric breath- 
ing modes - often soften across electronic transitions. For 
this reason, monitoring the evolution in Raman modes 
associated with breathing-type octahedral distortions can 
shed light on the origin and onset of charge ordering 
metal-insulator transitions in bulk perovskites or super- 
lattices formed by interleaving such materials.^ In some 
orthorhombic structures, the Jahn- Teller activated modes 
are compatible with the crystal structure, i.e., there is 
no symmetry change, due to two different types of B— 
O bond lengths induced by the mode. In contrast, the 
identical B-0 bonds present in rhombohedral structures 
would require that the Jahn- Teller activated modes reduce 
the crystal symmetry — the appearance of these modes 
in a Raman spectrum would reveal a symmetry-breaking 
phase transition. 

The combined sensitivity of Raman spectroscopy to 
symmetry, structural lattice dynamics and correlated elec- 
tronic and polar degrees of freedom, and versatility (e.g., 
facile application of fields and temperature) makes the 
technique well-suited for probing the capabilities of newly 
identified octahedral engineering approaches!^ Reports 
on Raman scattering selection rules and experimental 
work on then, however, are rare. In this work we present 
a comprehensive list of the normal modes, their vector 
displacement patterns, and the associated selection rules 
for the 15 octahedral rotation patterns (Glazer systems). 
We have systematically generated 705 schematic atomic 
displacement patterns for the normal modes of all 15 (14 
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FIG. 1. (Color online) The crystal structure of perovskite 
oxides with the ABO3 stoichiometry (a) The same structure as 
in (a) but showing the BOe octahedral network (b). Schematic 
illustration of the the a°b + b + structure (c), indicating that 
adjacent octahedra within the same plane rotate in opposite 
sense to maintain corner connectivity. Schematic illustration 
of the (d) a°a°c + and (e) a°a°c~ rotation patterns. The in- 
phase rotations and out-of-phase rotations of the octahedra 
about the z-axis are discernible in (d) and (e), respectively. 
[Note that for clarity the A-site cations and the oxygen atoms 
are omitted from panels (c), (d) and (e)]. 



rotated + 1 unrotated) Glazer tilt systems. Furthermore, 
we have shown how some recent experimental findings can 
be analyzed using our results and the schematic diagrams 
to draw decisive conclusions. Our results should be appli- 
cable immediately to the complete range of bulk ABO3 
perovskites with octahedral tilting. We further anticipate 
that equipped with comprehensive knowledge of the par- 
ent materials one should be able to analyze and rationally 
design more complex structures (e.g., epitaxial films and 
superlattices) composed of these ABO3 perovskites. 



II. STRUCTURAL DISTORTIONS IN 
PEROVSKITES 

The ABO3 type perovskite is built from basic building 
units of BOq octahedra that are corner-connected with 
additional charge balancing A-site cations occupying inter- 
stitial positions (Fig. [I]). The ideal, aristotype structure, is 
simple cubic with space group Oi or equivalently Pm3m 
in the Herman-Maguin notation. The A and B cations 
occupy sites with full cubic symmetry (point group Oh), 
whereas the O anions occupy sites with symmetry. 
In the ideal structure the ionic radii rvt, and ro are 
geometrically related by a characteristic tolerance factor 
t = (r A + r )/V2(r B + r ) = l- 17 The majority of the 
perovskites, however, are not cubic. The ionic radii of 
A, B and O rarely satisfy the condition that t = 1. As 
a result, the ideal cubic structure undergoes distortions: 
Typical distortions are rotations of the B0 6 octahedra 
about one or more high-symmetry axes, (anti)-parallel 
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displacement of the A or B cations away from their cubic 
site symmetry, or some combination thereof. 

The ABO3 type perovskites can also deviate from the 
ideal cubic structure by distortions to the rigid B0 6 oc- 
tahedra, i.e., by changes in B-O bonds lengths. These 
type of distortions are, however, mostly driven by va- 
lence fluctuations or orbital degeneracies, rather than 
steric constraints. Changes in B-0 bond lengths are 
possible and expected in systems with d 4 or d 7 electronic 
configurations on the B atoms, leading to Jahn- Teller 
antisymmetric breathing of the BO$ octahedra. Using 
group theoretical techniques, Stokes et al. found that 
among the 15 perovskite systems with simple rotation 
patterns, only one tilt pattern (space group P^/nmc, 
no. 137) requires distortions to the octahedra, i.e., it can- 
not be accommodated without changes in the B-0 bond 
lengths. 18 In general, however, departure from the cubic 
symmetry to a structure of lower symmetry will occur 
with both rigid rotations of octahedra and bond length 
distortions as the latter, secondary modes, often produce 
further energy stabilization of the crystal structure. 

A. Octahedral Rotation Syntax 

The classification of perovskite crystals according to the 
rotatio n patt erns of rigid BO§ octahedra is attributed to 
GlazerP^iSlEU The corner-connectivity of the BOq octa- 
hedra in the perovskite structure constrains the rotations 
of the adjacent octahedra. The rotations of an octahe- 
dron about a given axis requires that successive octahedra 
along directions perpendicular to that axis rotate in the 
opposite sense about that axis. Fig. file) shows that a 
positive rotation (e.g., counterclockwise) of an octahedron 
about the z-axis with magnitude a results in a negative 
(clockwise) rotation about z of equal magnitude of the 
adjacent octahedra located along the x and y directions. 
However, successive octahedra along the rotation axis, 
z, can have either the same or opposite rotation sense 
[Fig. gd) and 1(e)]. 

To simplify and generalize this description, Glazer in- 
troduced a notation whereby rotations in perovskites are 
classified according to how adjacent octahedra rotate 
along a particular Cartesian axis passing through the B- 
site: a~(a + ) indicates out-of-phase (in-phase) rotations 
while a denotes no rotations. For example, in a per- 
ovskite with the rotation pattern a + b~b~ the adjacent 
octahedra along the x-axis rotate in phase (e.g., all clock- 
wise) and the adjacent octahedra along the y and z axes 
rotate out of phase by the same angle b. Examples of the 
two single tilt systems, a°a°c + and a°a°c~, are shown 
schematically in Fig.[TJd) and 1(e). 

Rotations of the octahedra double the repeat distances 
perpendicular to the rotation axis and can result in unit 
cell vectors inclined to each other by angles other than 
90° due to ferroelastic strains.^ The transition into the 
distorted phase changes the number of formula units in 
the unit cell and the full space group symmetry of the 



crystal system. By considering all possible rotation senses 
and amplitudes along various Cartesian directions, Glazer 
determined that there are a total of 23 different possible 
rotation patterns. (They are also called tilt systems in 
the literature.) Some of the tilt systems yield identical 
space group symmetries: Group-theoretical analysis by 
Stokes and HatctP^ shows that of the 23 tilt systems, 
15 simple rotation patterns produce 15 centrosymmetric 
space groups. Table I shows these 15 tilt systems with the 
appropriate Glazer symbols, the space group notations 
(and numbers in parentheses) , the number of formula units 
(Z), the point symmetry, and known example compounds. 
The eight tilt systems that are absent from the group 
theoretical analysis were found to have higher symmetry 
than the corresponding space groups. Exhaustive studies 
by Woodward show that 12 of the 15 tilt systems reported 
by Howard and Stokes naturally exist (Table I). None of 
the eight Glazer tilt systems omitted by Howard and 
Stokes have been experimentally observedP 

III. CRYSTAL SYMMETRY: POINT GROUPS 
AND SPACE GROUPS 

Before describing the method for normal mode and 
selection rule determination in the Glazer systems, we de- 
scribe the general methods available to determine normal 
modes and selection rules in crystals. We begin by provid- 
ing a brief discussion of point groups, space groups and 
the relevant symmetry issues pertinent to this analysis. 

Bravais Lattice. — The problem of classifying crystal 
structures starts with defining the Bravais lattice, which 
is an infinite array of discrete points that appears iden- 
tical from every point of the array. The description of a 
physical crystal requires describing the underlying Bravais 
lattice and the arrangement of atoms or molecular build- 
ing blocks within a primitive cell. Each Bravais lattice 
is characterized by a set of rigid operations that takes 
the lattice into itself. This set of operations is known as 
the symmetry group of the Bravais lattice and includes 
all translations through the lattice vectors and the point 
symmetry operations, which are rotations, reflections, 
and inversions, or some combinations of these, described 
below. 

The point symmetry operations form a subset of the full 
symmetry group and are sufficient to construct the seven 
crystal systems — cubic, hexagonal, tetragonal, orthorhom- 
bic, rhombohedral, monoclinic and triclinic. If, however, 
the full symmetry group, which includes translations, 
are included, some of the crystal systems are expanded 
and fourteen distinct groups result — the fourteen Bravais 
lattices. For example, the cubic crystal system, after 
considering translations, expands to include simple cubic 
(sc), body centered cubic (bcc), and face centered cubic 
(fee). 

Point Groups. — A physical crystal is obtained when 
an object (e.g., a molecule) is attached to the Bravais 
lattice. Although in principle objects can have many 
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TABLE I. The Glazer systems found in perovskites. The space group symmetries, point symmetries, formula units, and 
experimental compounds, where available, are enumerated for each rotation pattern. The references for the experimental 
compounds can be found in reference 3 and are not be repeated here."-" indicates no known example. 



Tilt 


Space 


Point 


Known 


Systems 


Group 


Symmetry 


Example 


1. a°a°a 


Pmlm (221) Z = 1 


m3m (Oh) 


SrTiOs 


2. a°a°c+ 


P4/mfem (127) Z = 2 


4/mmm {Dah) 


CsSnI 3 
(351 - 425 K) 


3. a°b+b+ 


I4/mmm (139) Z = 8 


4/mmrn (Dih) 


_ 


4. a + a + a + 


7m3 (204) Z = 8 


m3 (T h ) 


Cao.25Cuo.75Mn03 


5. a + b+c+ 


I/mmm (71) Z = 8 


mmm (£>2h) 




6. a°a°c- 


/4/mcm (140) Z = 4 


4/mmm (D^h) 


SrTi0 3 
(<110 K) 


7. aPb'b- 


/mma (74) Z = 4 


mmm (D2h) 


PrAlOa 
( 151 - 205 K) 


8. a~a~a~ 


P3c (167) Z = 4 


3m (I>3h) 


LaNiOs 


9. a°b~c- 


C2/m (12) Z = 4 


2/m (C 2h ) 


PrAlOs 
(<135 K) 


10. a-ft- fe- 


C2/c (15) Z = 4 


2/m (C 2h ) 


- 


ll. a~b~c~ 


PI (2) Z = 2 


i (CO 


W0 3 
( 230 - 300 K) 


12. a°b+c- 


Cmcm (63) Z = 8 


mmm (D21) 


SrZrAlOs 
( 973 - 1103 K) 


13. a~b~c~ 


Pnma (62) Z = 4 


mmm (£>2h) 


LaMnG-3 


14. a + b~c- 


P2i/m (11) Z = 4 


2/m (C 2/t ) 


GaLiBr 3 


15. a + a + c~ 


P4 2 /nmc (137) Z = 8 


4/mmm (-D 2 /t) 


CaFeTi 2 3 



configurations and infinite number of point symmetry 
groups, the crystalline environment puts severe limitations 
on the allowed configurations. From the infinite number 
of three dimensional point groups, only those where the 
rotational symmetry C n has n = 1, 2, 3, 4, and 6 are 
possible in a crystal. The crystallographic point groups 
are constructed out the seven crystal systems (Bravais 
lattice points groups) by systematically considering all 
possible ways of reducing the symmetry of the objects 
characterized by these groups. The conditions yield the 
set of 32 crystallographic point groups. For example, the 
cubic system is expanded to five groups which are, in 
Schoenflies notations, Oh, O, Th, Td and T. 

Space Groups. — When a crystal is formed by inserting 
a motif into a Bravais lattice, the full symmetry of the 
crystal depends on both the symmetry of that object and 
that of the Bravais lattice; the full symmetry group is 
called the space group of the crystal. If the object has 
spherical symmetry, the total number of space groups 
would be fourteen, i.e., the same as the number of Bra- 
vais lattices. However, objects rarely have the spherical 
symmetry, making the number of possible crystal space 
groups significantly larger. Some of the space groups are 
constructed by placing an object with the symmetries of 
each of the thirty-two point groups into those fourteen 
Bravais lattices within the same groups. For example, 



as shown earlier, there are three Bravais lattices of the 
cubic system and there are five crystallographic points 
groups belonging to the cubic system. To first order then, 
there are a total of 15 space groups belonging to the cubic 
system. In this way, one arrives at a total of 73 space 
groups, known as the symmorphic space groups. 

There are additional considerations, however, in the 
case of the space groups that do not result from simple 
compounding of the Bravais lattice systems with the 
respective point groups. In the space group, a point 
is considered invariant if it is either left in place by a 
particular symmetry operation or is carried over to a 
position in an adjacent unit cell, which it could have 
reached by a simple translation of one unit cell. Also, the 
presence of different atomic motifs within the crystal and 
their different sets of symmetry operations means that 
the rotational axes and the symmetry planes need not 
coincide at a common point. As a result more than one 
set of operations may fall within the same class. 

For these reasons, two additional types of compound 
symmetry operations appear in space groups: the screw 
axis and the glide planes. A screw axis operation produces 
a rotation followed by a translation along the axial direc- 
tion while a glide plane operation yields reflection across 
a plane followed by translation along that plane. With 
the inclusion of the screw axis and the glide plane, 157 
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new space groups, known as the non-symmorphic group, 
emerge, yielding a total number of 230 space groups. The 
Glazer systems for perovskites with octahedral rotations 
are found in 15 (13 symmorphic and 2 non-symmorphic) 
of the general 230 available space groups. 

Site Symmetry. — Within each unit cell defined by a 
crystallographic space group, there are special sets of 
positions in which each point of the set has identical 
surroundings. These points, or sites, and the set of sym- 
metry operations that leave them invariant define the 
site group of the set of points; one point belonging to 
this set is called a Wyckoff position. The site group is a 
subgroup of the space group and is isomorphous with one 
of the 32 crystallographic point groups compatible with 
the space group. The number of equivalent points n in a 
site group is equal to the order of the factor group (space 
group) H divided by the order of the site group h, i.e., 
n = H/h. For example, in the simplest Glazer tilt system 
without any octahedral rotations, a°a°a° (space group 0\ 
or Pm3m), the A-site atom belongs to the site group la, 
with site symmetry Oh, the B-atom belongs to the site 
group lb with site symmetry Oh, and the three oxygen 
atoms belong to the site group 3c with site symmetry 
D^h The importance of site symmetry becomes evident 
in the next section where we discuss the various methods 
of normal modes determination in a crystal. 

Caveats. — A relevant issue in the discussion of the space 
groups and its relevance to normal mode enumeration, 
which is the aim of this work, is the relationship of the 
primitive cell of a crystal to that of the crystallographic 
cell, as that interrelationship is not always straightforward. 
A primitive cell is the smallest volume from which the 
entire crystal can be reproduced by simple translations 
through the primitive vectors. It contains a single lattice 
point which is invariant under all rotational operations of 
the space group. 

The normal mode determination for space groups with 
primitive lattices (designated by P in the space group 
label) is simple because in these cases the crystallographic 
cell corresponds to the primitive cell. In the case of space 
groups with non-primitive (centered) lattices, the usual 
construction of primitive cell is less useful because it does 
not easily reveal the full symmetry of the structure. For 
centered lattices, the crystallographic cell is constructed 
from an integral number of the primitive cells obtained by 
translating the primitive cell lattice point in a prescribed 
manner to reveal the full symmetry of the structure. The 
translation vectors are then redefined to correspond to 
the expanded, non-primitive, cell. Although the primitive 
cell is rarely used to describe the crystal structure, the 
normal mode evaluation relies on the primitive lattice and 
not the most convenient setting. Care should therefore be 
exercised when determining the character tables and the 
normal modes of space groups with non-primitive lattices. 
We illustrate these issues later as some of the 15 Glazer 
tilt systems, are described by body- and base-centered 
cells. 



IV. NORMAL MODE DETERMINATION 

A. General Approaches 

Several methods have been developed to analyze the 
primitive unit cell and determine the selection rules for 
the first order phonon spectrum, i.e., one phonon for 
one photon under non-resonance conditions, at the zone 
center (r-point). We describe those approaches here 
briefly, with special emphasis on the Nuclear Site-group 
Analysis (NSA) method, which we use to enumerate the 
normal modes for perovskites crystals with octahedral 
rotations. In the next section, we discuss the application 
of the NSA process in this context in detail. 

(i) Factor group analysis: In this method the symmetry 
properties of the crystal are determined by studying 
the effect of each symmetry operation in the factor 
group (s pace group) on each type of atom in the 
unit cell! 24 * 25 ! DeAngelis et al. have explained this 
method with extreme care in Ref. 1231 

(ii) Molecular site group analysis: In this approach the 
symmetry properties are determined first for the 
asymmetric unit of interest (e.g., a molecule) as an 
isolated speciesP^ These results are then analyzed 
in terms of site symmetry and finally in terms the 
factor group (space group) symmetry. 

(Hi) Nuclear site group analysis (NSA): In this tech- 
nique a site-symmetry analysis is carried out on 
every unique set of atoms in the unit cell. Tables 
are then constructed to allow symmetry and selec- 
tion rule determination based on the irreducible 
representations (irreps) of the space groups without 
the construction of a correlation diagram! 27 ! 28 ! 

The factor group analysis method is the most rigorous 
and requires a thorough understanding of the unit cell. 
Furthermore, since the symmetry analysis is performed on 
each individual nucleus, a separate symmetry generating 
coordinate (e.g., a Cartesian coordinate vector located on 
each nucleus) is needed for each nucleus. 

The two site group analysis methods are often used 
in studies on molecular crystals such as carbonates, sul- 
fates, phosphates and the like, where the forces between 
molecules are considerably weaker than those between 
atoms inside the molecule. However, the molecular site 
analysis method is a multistep process where the normal 
modes for the ion groups are determined first, followed 
by the rotational (libration) and the acoustic motions 
of the molecules. Since for most practical purposes the 
ABO3 perovskites can be considered as a flexible network 
of corner connected BOq octahedra, the site group anal- 
ysis methods, especially the NSA method is used more 
extensively. We follow this route here. Indeed, both the 
rigorous factor group analysis method of DeAngelis et al. 
in Ref. 25 and our nuclear site group analysis method 
(Table 13) produce normal mode decompositions with 
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identical numbers and symmetries for the Pnma space 
group. 

To understand the approach in the NSA method, con- 
sider for example, a set of atoms with site group g inside 
the unit cell of a crystal of space group G. If the order 
of G is H and that of g is h, then n = H/h, where n 
is the number of equivalent points (e.g., atoms) in the 
site group g. Since g is a subgroup of G, each irrep of 
g will map onto one or more of the irreps of G. In the 
NSA method only the vector irreps of the site group g are 
chosen and these are then mapped onto the irreps of the 
crystallographic point group of the space group G. The 
total number of normal modes of a crystal with a partic- 
ular space group can finally be obtained by considering 
all the site groups and occupied Wyckoff positions and 
taking their algebraic sum. 



B. Methodology for Perovskites with Octahedral 
Rotations 

We use the Nuclear Site Group Analysis (NSA), method 
(Hi), to enumerate the normal modes of ABO3 perovskite 
oxides with only distortions characterized by rigid rota- 
tions of BOq octahedra. We begin with the prototypi- 
cal aristotype cubic polymorph, e.g., room-temperature 
SrTiC>3, with the symmetry 0\ (PmZm). To generate 
the rotationally distorted phases, we then freeze-in each 
specific symmetry lowering rotation pattern (Table I, col- 
umn 1) into the cubic structure to obtain the space groups 
in column 2 of Table I. The rotations lead to a splitting 
of the oxygen Wyckoff orbits, which have been tabulated 
elsewhere!^ 

Knowledge of the occupied cation and anion Wyckoff 
positions in the standard notations is essential to the 
NSA method, because they are directly related to the 
site symmetry of the atomic species. The normal modes 
of each atom are identified by consulting the Interna- 
tional Tab le of X -Ray Crystallography (ITXRC) and/or 
Porto et al! 28 * 3 ° l for the character tables of the site groups 
and of the space group. These elements are required 
for determining the vector irreps of the site group and 
subsequently mapping these to the irreps of the space 
group. We repeated this analysis for all atomic sets until 
the total number of normal modes for the asymmetric 
unit was obtained. The data are presented in Tables III - 
XVII found in the appendix. 



C. Example: Normal Mode Determination of the 
a + b~b~ Rotation Pattern 

We now describe this process in detail for the or- 
thorhombic a + b~b~ rotation pattern, Glazer tilt system 
13 (Table III), which LaMn0 3 is known to exhibit.^ 
When the tilt system a + b~b~ is applied to the cubic 
structure the space group transforms from cubic Pm3m 
to orthorhombic Pnma. The displacement pattern of 



TABLE II. Occupied Wyckoff positions and site symmetries 
for Pnma orthorhombic (space group 62, a = 5.743, b = 7.695, 
c = 5.537 A) LaMnOs from Ref. [31] 





Wyckoff 


Site 




Position 




Atoms 


Orbit 


Symmetries 


X 


y 


z 


La 


4c 




0.55 


1 

4 


0.009 


Mn 


4a 


a 











0(1) 


4c 


cr 


-0.011 


1 

4 


-0.071 


0(2) 


8d 


Ci 


0.309 


0.039 


0.225 



oxygen atoms also leads to an approximate \[2 x \/2 x 2 
increase in the cell dimensions a, b and c, respectively, pro- 
ducing a unit cell with four ABO3 formula units (Z = 4), 
as enumerated in column 2 of Table I. In Schoenflies nota- 
tion, the Pnma space group is equivalent to D\^, making 
the relevant point symmetry mmm or D 2 h (Table 13). 
After obtaining the lower symmetry structure ( Table II I , 



the Wyckoff positions can be read from the ITXRC, or 
by using crystallographic vi sualiz ation software such as 
vesta and CrystalMakerP^ Note that in the Wyck- 
off orbit column, the numeral preceding the Wyckoff label 
indicates that all of those symmetry equivalent atoms 
occupy that particular Wyckoff position. Chemically 
equivalent atoms can occupy crystallographically unique 
Wyckoff positions in crystals. This becomes apparent in 
LaMnOa with the a + b~b~ rotation pattern: The oxygen 
atoms occupy sites 4c and 8d ( Table II ) , whereas in the 



ideal cubic case only one Wyckoff position is necessary to 
describe all oxygen atoms in the crystal structure. 

The eight oxygen atoms parallel to the ac plane oc- 
cupy Wyckoff positions d and possess site symmetry 
Ci, whereas the four oxygen atoms along the b axis oc- 
cupy Wyckoff positions c and possess site symmetry Cf z . 
This "splitting" of the Wyckoff orbits of the O atoms 
is a signature of the structural transition on going from 
Pmim^-Pnma, which can be exploited in Raman spec- 
troscopy studies to recognize the onset of the transition 
by monitoring the change in number of peaks (modes), 
intensities, etc., with external stimuli. 

Consider, the four oxygen atoms occupying the Wyckoff 
position 4c and site symmetry C s , the irreps of w hich 
can be found in the ITXRC and/or Porto et al. } Mm The 
vector irreps of C s are A' and A" . When these are mapped 
on to the irreps of the crystallographic point group D2/1 of 
Pnma, we obtain the contributions of these four oxygen 
atoms to the normal modes of LaMnOa . This is repeated 
for all the atomic species in the unit cell. 

The first column of Table XV shows all available normal 
modes for the tilt system a + b~b~ in the absence of any 
other atomic distortions (cation displacements or Jahn- 
Teller distortions of the ligands). The values in the main 
body of the table shows the number of symmetry modes 
in which the atoms participate. This type of presentation, 
instead of simply writing down all the modes, provides 
physical meaning to the modes. As shown in Table XV, 
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the 4 La atoms have site-symmetry C'j. z and contribute 
the following modes: 

2A g © A u © 2Bi 9 © B lu © 2.82;, © £ 2u © B 3g © 25 3 „ . 

The 4 Mn atoms with site-symmetry Ci contribute: 

3A U © 3B lu © 35 2u © 3B 3u , 

while the 4 O(l) atoms with site-symmetry C* z con- 
tribute: 

2A g ®A U ® 2B lg © B lu © 2^23 © £ 2u © B 3g © 25 3u , 
and the 8 0(2) atoms with site symmetry C\ contribute: 
3A g © 3A U © 3Bi g © 3B lu © 3B 2s © 3B 2u © 3£ 3g © 3B 3u . 

The last two columns in Table XV provide the mode 
selection rules in standard notations. The acoustic modes 
belong to representations which contain the translation 
vectors, symbolized by T. After subtracting out the three 
acoustic modes available to three-dimensional crystals, the 
remaining modes labeled by T are IR active. The selection 
rules for the Raman modes are shown separately in the 
last column with the polarizability tensor components, 
denoted by and the appropriate Cartesian subscripts. 
OLij are also known as susceptibility derivative tensors in 
the literature. 

For Puma LaMn0 3 , Table XV reveals that one of each 
of the Bi ui B 2u and B 3u modes account for the acoustic 
(translational) modes of the crystal. The remaining 9Bi u , 
3i?2u and 8B 3u modes are IR-active, and the 7A g , 5B lg , 
7-B 2s , and 5£? 3ff modes are Raman active. The normal 
modes that are neither IR active, nor Raman active are 
considered silent. In the present case, the silent modes 
transforms with A u symmetry. 

For completeness, we also provide the rotation (axial) 
vectors compatible with the appropriate modes. For ex- 
ample, as shown in Table XV, the rotation R z has the 
same symmetry properties as that of the A g modes. In 
general, the rotations have the same symmetry proper- 
ties as those of the difference between anti-symmetric 
Raman tensors: R z transforms as a xy — a yx . For Ra- 
man scattering that uses excitation frequency far from 
electronic resonance, these differences are zero since the 
polarizability tensors are symmetric. In resonance Raman 
experiments, however, asymmetries in the polarizability 
tensors are introduced and then such modes, i.e., modes 
with R x , R y , R z labels could become Raman active. 

The last row of each table shows the total number of 
modes for each rotation pattern. Note that modes with A 
and B spectroscopic labels are non-degenerate, whereas 
modes with E (T) labels are doubly (triply) degenerate. 

Finally, to complement these tables we have developed 
a graphical description of each mode by using the charac- 
ter projection method. The complete compilation of all 
705 of these diagrams is available by accessing: |LINK A 



few of the representative diagrams are shown in |Figure 2 



The vector displacements without such diagrams, corre- 
sponding to all normal modes for each Wyckoff position 
available to the 230 space groups, can also be obtained us- 
ing the Bilbao Crystallographic Server and the web-based 
program SAM 



V. DISCUSSIONS AND APPLICATIONS 



A. Raman Mode Bookkeeping: How many modes 
are there? 



One way to check the accuracy of the normal mode 
tables we have generated for the ABO3 perovskites with 
octahedral rotations (Tables III-XVII) is to count the 
number of normal modes enumerated for each crystal 
system. If there are p > 2 different atoms in a three- 
dimensional primitive cell, there should be 3p branches 
to the phonon dispersion relation, and consequently, 3p 
modes at the T point. Of these 3p modes, 3 are acoustic 
(one longitudinal and two transverse) modes and the 
remaining, 3p — 3, are optical modes. For space groups 
with a primitive lattice, as is the case for the Glazer tilt 
systems 1, 2, 8, 11, 13, 14, and 15, the Raman mode 
enumeration is straightforward. 

For the Glazer systems with non-primitive lattice, care 
must be taken to avoid redundancy in the enumeration of 
the normal modes. As discussed earlier, for these systems 
the usual construction of the primitive cell is less useful 
because it doesn't usually reveal the full symmetry of 
the structure. In these cases the crystallographic cell is 
constructed by tiling an integral number of primitive cells 
through lattice translations of the primitive cell lattice 
point. Since the translation group transforms as the 
identity element of the space group, the distribution of 
degrees of freedom among the space group irreps does 
not change. However, this introduces a redundancy in 
the value of the characters for each symmetry: Each set 
of equipoints in the primitive cell is duplicated for the 
additional lattice points in the expanded crystallographic 
cell. Since redefining the unit cell vectors does not affect 
the rotational symmetry, the redundancy is uniformly 
distributed over the characters of all symmetry classes of 
the space group. 

To recover the primitive (p) cell characters, one should 
divide the conventional (c) cell characters by the lat- 
tice point multiplicity, i.e., \p — Xc/£, where £ is the 
number of lattice points (or primitive cells) in the full 
crystallographic cell: £ = 4 for face-centered cell, £ = 3 
for hexagonal or rhombohedral cells and £ = 2 for both 
body centered cells and base centered cells.^ This con- 
sequently yields the number of normal modes being an 
integral fraction of 3q, where q is the number of atoms in 
the conventional cell. The Glazer systems 3, 4, 5, 6, and 
7 are body centered and as such the number of modes 
is 3q/2. The Glazer system 8 is rhombohedral and as 
such the number of modes is 3q/3. The Glazer systems 9, 
10 and 12 are base centered and as such the number of 




FIG. 2. (Color online) Schematic atomic vibrational patterns of selected normal modes for a few octahedral rotation patterns. 
The symmetry mode symbols (irreducible representations) and the atoms participating in the normal modes are also provided. 
The complete database for all normal modes in each Glazer tilt system can be viewed online by accessing: |INSERT TARGET| 
IADDR.ESSI 



modes is Zq/2. 



B. Applications: Phase Transitions in Bulk 
Perovskites 

Having ensured that the number of modes for each 
system is accurate, we now provide some illustrative ap- 
plications and discuss a few limitations of the constructed 
tables and the schematics. Here we show how the Raman 
normal modes determination of the various octahedral 
rotation patterns in perovskites provides a manner to 
explore antiferrodistortive phase transitions. 

CaTiC>3 is a widely used electroceramic and has been 
recognized as a suitable material for the immobilization 
of high-level radioactive wastesP^ Consequently, the prop- 
erties of CaTi03, especially its crystal structures over 
different temperature regimes, have lately been the fo- 
cus of intense studies.^ It undergoes a series of struc- 



tural transitions driven by the octahedral rotations (Fig 
ure 3 1 Under ambient conditions, CaTi03 is orthorhom 
bic (space group Pnma) and exhibits the a + b~b~ rota- 
tion pattern. Upon heating, the structure transforms to 
tetragonal (space group I4/mcm, Glazer system a°a°c~) 
at 1520 K. Note that some studies suggest that between 
the tetragonal and the orthorhombic phase there is an 
additional structure with orthorhombic symmetry, Cmcm, 
that exhibits one in-phase and one out-of-phase rotation 
(a°b + c~). Finally near 1645 K, CaTi03 transforms into 
the ideal cubic PmZm structure without any octahedral 
rotations (a°a°a°). 

Unpolarized Raman studies enable identification of 
the nature and onset of the phase transition from the 



' 1 




a+b-b- 




(Pnma) 








297 K N 


V 



1520 K 



(1 4/ mem) 



1645 K 



(Pm3m) 




FIG. 3. Sequence of structural phase transitions and the 
change in octahedral rotation patterns that occurs with in- 
creasing temperature in perovskite CaTi03. The orthorhombic 
Cmcm phase (linked by broken lines) appears (or is sometimes 
absent) experimentally as an intermediate phase between 
Pnma and 74/mcm at intermediate temperatures. Transi- 
tions required by Landau theory to be first (second) order are 
distinguished by heavy (light) lines. 



octahedrally rotated structures to the cubic structure 
upon heating. The relevant information for the analysis 
of the rotation-driven structural transitions 



~b~b~ 



a a°c 



->■ a°a°a 



are found in Tables XV, XIV, VIII and III, respectively. 
First, comparing these tables reveals that all of the rota- 
tionally distorted perovskite polymorphs have first-order 
Raman active vibrational modes, whereas all zone-center 
modes in cubic (aPa°a°) structure are Raman inactive. 
Accordingly, unpolarized Raman spectra of the tilted 
structures should exhibit first-order Raman features, pos- 
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sibly superimposed on broad second order features, where 
as the cubic should not. (Note that in practice, how- 
ever, at the transition to the cubic structure the first 
order Raman features could give way to broad second 
order features — such effects are seen in prototypical cubic 
perovskite SrTi0 3 P) 

Polarized Raman studies can be used to determine the 
phase transitions among the tilted structures. Raman 
scattering cross sections {<jr) for Stokes scattering is re- 
lated to the polarizations of the incident (e\) and scattered 
(e J s ) radiations and the susceptibility (derivative) tensors 
by the following formula: 

a R oc l&gre^ij'l 2 i 

where onj are the components of the susceptibility (deriva- 
tive) tensors appropriate for a given mode, shown in the 
last column of the Tables III-XVII. Note that summation 
over the repeated cartesian indices (i, j) is implicit in the 
above equation. 

As an application of this, let's consider Raman spec- 
troscopy in the back scattering geometry for CaTiC>3. 
The standard scattering geometry for this experiment is 
z(ij)z, which is shorthand for ki(e},e J s )ks, where ki and 
ks are the wave vectors of the incident and scattered 
light, respectively, and e} and e^ are their corresponding 
polarization vectors. Since both k$ and ki are in the z 
direction, the polarization vectors of both the incident 
and scattered radiation are in the xy-plane, resulting in 
Raman cross section of: 

a R cc \e x s e x a xx + e%e v I a xy + e y s e x a yx + e y s e y jayy\ 2 . 

Consider now the transition from a + b~b~ — > a°a°c~. 
As seen from Tables XV and VIII, there are significant 
differences between the mode distributions in the two 
structures. In the a + b~b~ (Prima) structure, the Raman 
active modes are: 

7A g © 5B lg ffi 7B 2g © 5B 3g , 

whereas in the aPaPc~ (1 4 / 'mem) system, the Raman 
active modes are: 

Ag ffi Bg ffi 2^2,, © Eg ffi E U . 

Based on the larger number of Raman modes in Prima (24 
compared to the 8 in 14 /mem CaTiOs), there should be 
a readily detectable difference in the Raman spectra, i.e., 
one should expect significantly more first-order Raman 
features in the Prima structure than in 14/ mem, even 
in unpolarized Raman spectroscopy. The two phases can 
further be distinguished by performing a careful polarized 
Raman experiment, whereby an appropriate selection of 
the incident and scattered light orientation will enhance or 
suppress a first-order Raman mode. Performing Raman 
spectroscopy in the z(xx)z geometry will results in Raman 
cross sections of 

a R (x \e x s e x a xx \ 2 . 



This will only show modes of A g symmetry of which 
are seven in the Prima and only one in the 14/ mem 
structure. Similar results can be expected in the z(yy)z 
geometry. Further careful studies with the z(yz)y config- 
uration would show first order Raman features from the 
singly degenerate B 3g mode in Prima CaTi03, whereas 
an Eg doublet would appear in the 14/mcm structure. 

The controversial phase transition from Prima to 
Cmcm would be more difficult to determine by exam- 
ining the Raman spectra alone, since the Raman active 
modes present in the two structures are similar (cf. Tables 
XV and XIV) . To remedy the ambiguity, one could rather 
use the octahedral rotation angle, specifically the com- 
ponent about the x-axis, as an order parameter for the 
Prima — ¥ Cmcm transition, because it goes from a finite 
to null value as the rotation pattern changes, i.e. a + b~b~ 
to a°b + c~. Recently, Dubroka et al. and Abrashev et 
al. have shown how to correlate the Raman A g mode 
frequency with the tilt angle, and su bsequ ently that tilt 
angle with the free oxygen positions! 13 * 15 ! Their results, 
along with our schematic diagrams for the normal modes 
of the rotationally distorted perovskites and lattice dy- 
namic calculations could be used to asc ertain the phase 
transition between these two structures. ) 35 l 36 l 



C. Applications: Structure— property Relations in 
Perovskite Thin Films and Superlattices 

We expect this methodology to have higher impact in 
ABO3 perovskite systems realized via octahedral engineer- 
ing route whereby efforts are being made to exploit epitax- 
ial strains, interfacial octahedral coupling, and superlat- 
tices to stabilize structures inaccessible via conventional 
solid-state chemistry techniques. Recent computational 
studies predict selective control of the octahedra about 
different crystal directions by using epitaxial strain to con- 
trol the rotation-lattice coupling in epitaxial films |2MU 
Strain can therefore be used to directly tailor the flavor - 
magnitude and phase - of octahedral rotation patterns in 
thin films. 

The normal mode determination presented here, could 
be used to design and analyze ABO3 films with octahedral 
rotations. Consider for example, Lai/ 3 Sr 2 /3Fe03 (LSFO), 
which is rhombohedral (Glazer tilt system a~a~a~) in 
the bulk. Its normal modes appear in Table X. However, 
if an epitaxial film of LSFO is grown on cubic SrTi03, 
it would likely transform to a monoclinic phase with the 
a~b~b~ octahedral rotation pattern due to lattice mis- 
match between the film and the substrate in the epitaxial 
plane. Reduction of symmetry will lead to a splitting 
of some degenerate modes of the rhombohedral system, 
which can be determined by using the standard reduction 
formulae^! and our Table XII, i.e., the relevant normal 
modes for the monoclinic LSFO. By comparing the Ra- 
man spectra of bulk LSFO and the epitaxial film, it should 
be possible to evaluate the change in crystal structure for 
the rhombohedral to monoclinic transition. 
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Formation of a coherent perovskite-perovskite het- 
crointerface produces a geometrical constraint of the oc- 
tahedral connectivity across the interface in synthetic 
structures.^ The octahedral rotations present in a sub- 
strate, for example, can propagate into the near-interface 
region of a thin film: This was recently observed ex- 
perimentally using high-resolution x-ray diffraction and 
synchrotron radiation, in (LaNi03) n /(SrMn03)2 superlat- 
ticesP^ This so-called octahedral proximity effect remains 
to be harnessed for materials design, yet has the abil- 
ity to mediate electron-lattice coupling across interfaces 
and produce remarkable changes in the electronic prop- 
erties of ultrathin films and superlattices. Segal et al. 
have shown that the propagation of lattice vibrational 
modes (phonons) associated with the dynamic octahedral 
rotations causes drastic changes in the resistivity of a 
Lao.53Sr .47Mn03 (LSMO) films as the system is cooled 
through the a°a°a° — > a°a°c~ phase transition of the 
SrTiOa substrate. The proposed microscopic mechanism 
was attributed to evanescent cross-interface coupling be- 
tween the charge carriers in the film and the soft phonon 
mode responsible for the out-of-phase rotation of TiOg in 
SrTi03. In other words, enhanced electron-lattice cou- 
pling mediated through the correlated oxygen octahedral 
motions P-^ 

Other external stimuli besides temperature can pro- 
mote novel behavior mediated by octahedral rotations. 
For example, Caviglia et al. have shown that an opti- 
cal excitation in resonance with a stretching mode of a 
perovskite substrate can trigger a dynamic (transient) in- 
sulator to metal transition in NdNiC>3 epitaxial films below 
the bulk Neel temperature.^ Recent Raman spectroscopy 
studies on superlattices of superconducting YBa2Cu30y 
and magnetoresistive La 2 /3Sr 1 / 3 Mn03 below the critical 
temperature of the cuprate, reveal that the supercon- 
ducting state can renormalize the lineshape of the MnC>6 
octahedral rotation modes in the manganite, providing an 
unexplored platform to tailor electron-lattice interactions 
in correlated materials.^ 

The unique effects described here are all provided by the 
normal mode distributions within the materials and their 
coupling to the low-energy electronic structure. They 
could possess tremendous technological potential, yet the 
microscopic mechanisms in nearly all cases remain in- 
adequately understood and characterized. Much of the 
current results of computational studies, for example, 
await experimental verifications. Segal et al. suggested 
that the change in conductivity in their LSMO epitax- 
ial film could be due to the alteration of the B-O-B 
bond angles in the film by the soft phonon modes of the 
STO substrate — direct measurement of these variations 
remains to be seen. Note however that they also offered 
an equally likely scenario: The substrate could modify the 
B-0 bond lengths in the film, resulting in similar modifi- 
cations to the films conductivity, via a bond stretching 
mode. In the experiments by Caviglia et al., the authors 
acknowledge that the precise nature of the interaction 
between the LAO substrate and the NdNi03 film remains 



poorly understood and merits systematic studies. For 
that reason, the group offered an alternative explanation 
for the dynamic insulator-metal transition: The stretch- 
ing mode of the substrate could couple to a rotational 
mode (soft mode) of the film, which would also alter the 
films conductivity through a modification of the Ni 3d - 
O 2p angular orbital overlap. 

For such complex systems, our results presented here, 
especially the normal mode schematics, could be useful 
in drawing more decisive conclusions. By superimposing 
a unit cell of an epitaxial material on a unit cell of the 
intended substrate and inducing a particular normal mode 
in the substrate, once could examine geometrically and 
on symmetry grounds the extent to which the substrate 
modes couple to the film and induce changes in both the 
topology of the film's perovskite lattice and properties. 
Since our paper presents a complete list of normal modes, 
selection rules, and schematic diagrams of the bulk ABO3 
perovskites with octahedral rotations, careful implementa- 
tion of our results could guide research avenues in search 
of quantum phase transitions in artificial materials using 
optical excitation of lattice modes. 

D. Limitations 

We believe we would be remiss if we omitted a discussion 
of some limitations of the current paper and the methods 
used to derive the normal mode distributions. The normal 
modes in our paper were derived for simple, single phase, 
ABO3 perovskites with octahedral rotations. In mixed 
A- and J3-site systems, e.g. Lai-ySryMni^^A/^FeOs (M 
= Cr, Co, Cu, Zn, Sc, or Ga), the normal modes may 
be drastically different from their simple ^4-603 counter- 
part.^ In these cases the A and B site disorder suppresses 
the translational invariance that is a requirement for the 
range of applicability of the symmetry analysis. The re- 
sult being that the lineshapes of the modes with atomic 
displacement patterns from those sites could be broad- 
ened. Similar broadening could occur even in single phase 
perovskites grown in thin film form, where unintentional 
disorder and variability in the site occupation is induced 
preferentially during growth along various directions, i.e., 
in the plane of or perpendicular to the substrate. Such 
broadening can further lead to linewidths in which dis- 
tinct modes of different, but closely spaced, energies are 
indistinguishable. 

Second, our results were derived for bulk perovskites. 
Their direct application to epitaxial films and superlat- 
tices may not be immediately transferable. As mentioned 
earlier, the Glazer octahedral tilt patterns are derived 
by considering the octahedra as rigid units. However, in 
epitaxial films and superlattices, the octahedra are nec- 
essarily distorted. A consequence of such geometrically 
required distortions is that they can produce crystal struc- 
tures with space groups different than those analyzed here, 
i.e., those obtained from Glazer's analysis of possible com- 
binations of in- and out-of-phase rotations. Nevertheless, 
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our schemes would be a good starting point in analyzing 
these systems and deviations from ideality. For example, 
a comparison based on the Raman mode predictions pro- 
duced with our tables to the actual experimental results 
could help in understanding the perturbative nature of 
octahedral engineering in artificial perovskites. 

Lastly, our treatment omits the consequences of mag- 
netic ordering on the crystal symmetry, which is expected 
to become critical in perovskite superlattices formed with 
transition metal B-sites with open d-shell configurations. 
Several ABO3 type perovskites undergo paramagnetic 
metal to antiferromagnetic insulator transitions followed 
by charge ordering as they are cooled through their Neel 
temperature. The crystal symmetry operations in the 
antiferromagnetic state must leave invariant not only the 
positions of the ions but also their magnetic moments. 
For example, magnetic ordering in the rutile structure re- 
duces the space group symmetry from that of rutile with 
non-magnetic cations (P^jmnm) to Pnnm or D^j^ 
which in turn alters the Raman mode distribution. We 
hope our paper spawns additional developments in the 
experimental Raman studies of artificial perovskite ox- 
ides, theoretical approaches to predict the normal mode 
distribution in such materials, and means to overcome 
the current limitations discussed. 



VI. CONCLUSIONS 

Using the nuclear site group analysis method we have 
determined the complete set of normal modes and the 
associated selection rules for all 15 of the per- 
ovskites systems with octahedral rotations (the Glazer 
systems). For each of mode, we have produced a corre- 
sponding schematic diagram showing the vector displace- 
ment pattern of the atoms participating in the particular 
mode. The results of this analysis is a compendium of 
705 schematic diagrams, which are now web accessible 
through an online portal. 

We have shown how some recent experimental findings 
can be analyzed using our tables and the schematic di- 
agrams to draw decisive conclusions. We expect these 
tables and the schematic diagrams to be useful tools in 
analyzing the phase transitions, electron-lattice coupling 
and elucidating structure-property relationships in ABO3 
perovskites with octahedral rotations. Equipped with this 
database for the bulk perovskites, we suggest that one 
would be able to both analyze and rationally design func- 
tional artificial materials built out of these perovskite 
blocks and rotated octahedral units. 
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TABLE III. System 1— Normal Modes for the Glazer system a°b°c°. 



a°b°c 


A 
Oh{l) 




B 
Oh{l) 




D 4h (3) 


Selection Rules 


Tin 

Tin 


1 




1 




2 
1 


T 


Sum 


3 




3 




9 


15 






TABLE IV. System S 


! — Normal Modes for the Glazer system a°b°c + . 


a°b°c + 


A 

D' 2h (2) 


B 

C Ah (2) 




O(l) 
C ih (2) 


0{2) 
CS„(4) 


Selection Rules 


Alg 

Am 

A 2 g 

A 2u 

Big 

Dlu 
B 2 g 

B 2u 

E 3 

E u 


1 
1 

2 


1 
1 

2 




1 
1 

2 


1 

1 
1 
1 
1 
1 

1 

2 


Ot-Xy 

i^Rx^Ry) (Xxz)CKyz 

(T X ,Ty) 


Sum 


6 


6 




6 


12 


30 






TABLE V. 


System 3 


— Normal Modes for the Glazer system a°b + b + . 


a°b+b+ 


A(l) A(2) A{3) 
D 2h (4) T» 4h (2) D 4h (2) 


B 

Cm. (8) 


O(l) 
C?(16) 


0(2) Selection Rules 

C*,(8) 


Alg 
Alu 
A 2 g 

A 2u 

Big 
o 

Dlu 
B 2 g 

B 2u 

E 9 

E u 


1 

1 

2 


1 
1 


1 
1 


1 

2 
2 
1 
3 


2 
1 
1 
2 
2 
1 
1 
2 
3 
3 


1 i?, 
1 T z 

1 ^XX Oiyy 

1 

1 (R x ,R y ) (a xz ,a yz ) 

2 (T x ,T y ) 


Sum 


6 


3 


3 


12 


24 


12 60 






TABLE VI. System 4 


— Normal Modes for the Glazer system a + a + a + . 




A(l) 


A(2) 


B 


O 




Selection Rules 


a + a + a + 


Th(2) 


D 2h {6) 


S 6 (8) 


C s (24) 






A g 
A u 
E 9 
E u 
T g 
T u 


1 


3 


1 
1 
3 


2 
1 
2 
1 

4 
5 


R 
T 


Q^xa; "1" "1" (%zz 
(ttxj: Qyy 2q?22, ^ScXxx V^^i/y) 
C^:c:e ; OLxz ; C^yz 


Sum 


3 


9 


12 


36 




60 
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TABLE VII. System 5 — Normal Modes for the Glazer system a + b + c + . 





A(l) 


A(2) 


A(3) 


A(4) 


B 


O(l) 


0{2) 


0(3) 




Selection Rules 


a + b+c+ 


D 2h (2) 


D 2h (2) 


D 2h (2) 


D 2h {2) 


a (8) 


C*f z (8) 


C?»(8) 


cr(8) 






A g 












2 


2 


2 




(%xx i Qyy s ^23 


A u 










3 


1 


1 


1 






Big 












1 


2 


1 




&xy 


Biu 


1 


1 


1 


1 


3 


2 


1 


2 


T z 




B 2 g 












1 


1 


2 


Ry 


&XZ 


B 2u 


1 


1 


1 


1 


3 


2 


2 


1 


Ty 




Bag 












2 


1 


1 


R x 


Otyz 


Bi u 


1 


1 


1 


1 


3 


1 


2 


2 


T x 




Sum 


3 


3 


3 


3 


12 


12 


12 


12 




60 



TABLE VIII. System 6— Normal Modes for the Glazer system a°a°c". 



a°a°c- 


A 
D' 2d {±) 


B 

C 4 h(4) 


O(l) 

C 2 „(8) 


0(2) 
D 4 (4) 




Selection Rules 




Alg 






1 






a 


XX ~\~ &XX j 


A lu 




1 












A 2 g 






1 


1 








A 2u 


1 


1 


1 


1 








Big 






1 








O-xx Q-yy 


Bi u 






1 










B 2 g 


1 




1 








Oixy 


B 2u 
















Eg 


1 




1 


1 


Ry) 






E„ 


1 


2 


2 


1 








Sum 


6 


6 


12 


6 




30 





TABLE IX. System 7 — Normal Modes for the Glazer system a°b b . 





A 


B 


O(l) 


0(2) 




Selection Rules 


a°b~b- 


C a \(4) 


C? h (4) 


CI (8) 


Cf„(4) 






Ag 


1 




1 


1 




Otxx i Q-yy i Ot-zz 


A u 




1 


1 








Big 






2 




Rz 


Q-xy 


Biu 


1 


2 


2 




T z 




B 2 g 


1 




1 




Ry 


&xz 


B 2u 


1 


2 


1 




T z 




Big 


1 




2 




Rx 


a yz 


Bi u 


1 


1 


2 




T x 




Sum 


6 


6 


12 


6 




30 
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TABLE X. System 8 — Normal Modes for the Glazer system a a a . 



a a a 


A 
D 3 (2) 


B 
5 6 (2) 


O 

C 2 (6) 




Selection Rules 








1 




&xx , &yy , &zz 


A lu 




1 


1 






A 2g 


1 




2 






A 2u 


1 


1 


2 


T z 




E 9 


1 




3 


(Rx, Ry) 


(c*xx Qyy, &xy), (&xz, Qyz) 


E u 


1 


2 


1 


T T 

± x, J-y 




Sum 


6 


6 


18 




30 



TABLE XL System 9 — Normal Modes for the Glazer system a b c 





Ai 


B 


O(l) 


0(2) 


0(3) 




Selection Rules 


a°b- C - 


Cs (4) 


CM) 


C 2 (4) 


C a (4) 


C 3 (4) 






A 9 


2 




1 


1 


2 




Otxxi Qyyi OLzzi OL X y 


A u 


1 


3 


1 


1 


1 


T z 




By 


1 




2 


2 


1 




Otxx ■> Qyy 


B u 


2 


3 


2 


2 


2 


± xi J -y 




Sum 


6 


6 


6 


6 


6 




30 



TABLE XII. System 10 — Normal Modes for the Glazer system abb. 





A 1 


B 


O(l) 


0(2) 




Selection Rules 


abb 


C a (4) 


Ci (4) 


Ci(8) 


C 2 (4) 






Ay 


1 




3 


1 


Rz 


Otxx ■> OLyy , OL Z Z j <^XT/ 


A u 


1 


3 


3 


1 


T z 




By 


2 




3 


2 


Rx, Ry 




B u 


2 


3 


3 


2 


T X , Ty 




Sum 


6 


6 


12 


6 




30 



TABLE XIII. System 11 — Normal Modes for the Glazer system a b c . 





A 


B(l) 


5(2) 


O(l) 


0(2) 


0(2) 




Selection Rules 


a b c 


Ci (2) 


Ci(l) 


Ci(l) 


d(2) 


Ci(2) 


Ci(2) 






Ay 


3 






3 


3 


3 


R 


a 


A u 


3 


3 


3 


3 


3 


3 


T 




Sum 


6 


3 


3 


6 


6 


6 




30 
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TABLE XIV. System 12— Normal Modes for the Glazer system a b + c~. 





A(l) 


A(2) 


B 


0(1) 


0(2) 


0(3) 




Selection Rules 


a°b + c- 




C&(4) 


a (8) 


cr(8) 




Cf(8) 






A g 


1 


1 




2 


2 


1 




&xx 5 &yyi Oizz 


A u 






3 


1 


1 


1 






Big 


1 


1 




1 


2 


2 




&-xy 


B\ u 


1 


1 


3 


2 


1 


2 


T z 




B 2 g 








1 


1 


2 


Ry 


&XZ 


B 2u 


1 


1 


3 


2 


2 


2 


Ty 




Bag 


1 


1 




2 


1 


1 


Rx 


OLyz 


Bi u 


1 


1 


3 


1 


2 


1 


T x 




Sum 


6 


6 


12 


12 


12 


12 




60 



TABLE XV. System 13— Normal Modes for the Glazer system a + b~b 





A 


B 


Oil) 


0(2) 




Selection Rules 


a+b-b- 


C?*(4) 




or (4) 


d(8) 






Ag 


2 




2 


3 




Otxx i Qyy i Otzz 


A u 


1 


3 


1 


3 






Big 


1 




1 


3 




OL X y 


Blu 


2 


3 


2 


3 






B 2 g 


2 




2 


3 


Ry 


a xz 


B 2u 


1 


3 


1 


3 


Ty 




B 3 g 


1 




1 


3 


Rx 


Oty Z 


B: iu 


2 


3 


2 


3 


T x 




Sum 


12 


12 


12 


24 




60 



TABLE XVI. System 14 — Normal Modes for the Glazer system a + b^c'. 





A(l) 


A(2) 


5(1) 


B{2) 


Oil) 


0(2) 


0(3) 


0(4) 


Selection Rules 


a + b c 


(7.(2) 


C s {2) 


a (2) 


a (2) 


Ci(4) 


C.(2) 


C.(2) 


Ci(4) 




Ag 


2 


2 






3 


2 


2 


3 


Rz Qxx > &yy > &zz > &xy 




1 


1 


3 


3 


3 


1 


1 


3 


T 
± z 


Bg 


1 


1 






3 


1 


1 


3 


Rx,Ry Otxx,Ctyy 


B u 


2 


2 


O 


3 


Q 
O 


2 


2 


O 


Tx, Ty 


Sum 


6 


6 


6 


6 


12 


6 


6 


12 


60 
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TABLE XVII. System 15 — Normal Modes for the Glazer system a + a + c . 





A{1) 


A(2) 


A(3) 


B 


O(l) 


0(2) 


0(3) 




Selection Rules 


a + a + c 


D 2d {2) 


D 2d (2) 


C;(8) 


C?(8) 


a (8) 


(8) 








1 








2 


2 


1 






A lu 








3 


1 


1 


1 






A 2g 










1 


1 


2 






A 2u 


i 


1 


1 


3 


2 


2 


2 






Big 


i 


1 


1 




2 


2 


2 




(X-xx Ot X y 


-Dlu 








3 
O 


1 


i 
i 


o 
z 






B 2 g 










1 


1 


1 






B 2u 


i 






3 


2 


2 


1 






Eg 


2 


1 


1 




3 


3 


3 




{cX-xz i Qyz) 


E u 


2 


1 


1 


6 


3 


3 


3 


(T X ,Ty) 




Sum 


12 


6 


6 


24 


24 


24 


24 




120 



